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Abstract 



X 

In this paper, eit is shown for a single- j shell that the intrinsic excitations 
induced by the quadrupole-quadrupole interaction can be viewed as that due 
to the inverted surface-delta interaction whose multipole-pairing component 
with the largest spin is the dominant term and is attractive. This allow us 
to trace the origin of the recently introduced optimum pair description of 
rotation and leads naturally to analytical expressions for the rotational ener- 
gies for the various bands in a single j-shell with the quadrupole-quadrupole 
interaction 
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One of the central aims of nuclear structure physics is to understand how strong 
collectivity is manifested in the spherical shell model. Since straightforward application of 
the shell model is not possible for strongly deformed systems, a great deal of effort in the 
past forty years was devoted to the development of suitable truncation schemes to confront 
this problem 1 , 2 . 

It is well known now that under a strong monopole pairing interaction plus a weak 
quadrupole-quadrupole (QQ) interaction, monopole and quadrupole pairs appear to be a 
good truncation scheme. However, when the QQ interaction becomes dominant, it appears 
that higher multipoles, especially the hexadecapole pairs, must be added 3 . 

On the other hand, it has been recently been shown that the eigenstates of the pure QQ 
interaction in a single- j shell can be well described in terms of a collection of optimum pairs, 
using a novel approach which was proposed to describe nuclear states within a spherical 
shell model with an optimum pair basis 4 . This paper aims to trace the origin of such a 
description of rotation in a single- j shell with the strong QQ interaction. The physics behind 
this approach can be easily understood as follows. The major part of the QQ interaction in 
the single j shell can be expressed as a constant plus an L 2 term. The remaining parts can be 
simulated by an inverted Surface-Delta-Interaction (SDI), which will favor higher spin pairs. 
The constant plus the L 2 terms can be exactly diagonalized. This will give the rotational 
spectrum. The remaining part (i.e. the inverted SDI) can be approximately diagonalized by 
using first the standard many-body Tamm-Dancoff approximation, followed by the optimum 
pair approach. A by-product of this study is the derivation of analytical expressions for the 
rotational energies for various bands, for a system of any number of pairs moving in a singlej 
shell, although the number dependence is not understood at this point. In hindsight, this 
approach is possible because the intrinsic excitations induced by the QQ-interaction comes 
essentially from the inverted SDI. 

We shall illustrate our method via a simple solvable model: like-particles moving in a 
single- j shell. This example is chosen because the QQ-interaction can simulate the rotational 
excitation structure 2 . For a system of n/2 or N pairs of nucleons in single- j shell, the most 
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general Hamiltonian is 



H = e^j + l (at x a,) + ]T v / 2ATTC a (A a x A a )° (1) 

A 

with 

A A = ^ (a] x a]) A (2) 

and 

C A =< jjA || V || jjA > (3) 

where aj (aj) are the creation (annihilation) operators with energy ej. 

Using the equation-of- motion technique 5 , we shall employ the many-pair creation oper- 
ator as the excitation operator and the Tamm-Dancoff vacuum as the "true" vacuum. To 
portray a collection of multipole pair excitations, the excitation operator is expressed as a 
simple tensor product: 

fit = A t 1 Aj....A] v (4) 

where Aj stands for A^ , the tensor of rank L { . 

Following the standard process of linearization, the corresponding equation-of-motion be- 
comes 

[H, fit] = £ T fit (5) 

and the TDA energies are 

E t = J2NaC a (6) 

A 

Here N A is the number of such pairs with angular momentum A. 

It is known that the TDA is most valid when j is large and for interactions such as surface 
delta where the strengths of the multipole pairing components decrease monotonically. It is 
unsuitable for the identity interaction and the QQ force. For example, one can easily show 



that the TDA ground state energy for four particles in j = 21/2 orbit is only a third (a 
fifth) of the exact value for the identity (QQ) interaction. However, as will be shown later, 
the inadequacy of the TDA for these interactions can be salvaged by the introduction of an 
auxliary Hamiltonian to the QQ force which will render it as an " inverted" SDI interaction, 
one in which the strengths of the new multipole pairing components increase rather than 
decrase monotoniclly. It will restore the validity of the TDA. To this end, we shall introduce 
the following auxiliary Hamiltonian as 

H' = H QQ -aJ2^-lV (7) 

i<j 

where L is the angular momentum operator and 7^- stands for the identity interaction and 
L and are so chosen for their simplicity in representing an exact symmetry of the system. 
In fact, the auxiliary Hamiltonian of Eq.(7) has the same form as the original one given by 
Eq.(l) if the interaction matrix element, or multipole-pairing strength, are redefined as 

C' A = C K - a - 7 [A(A + 1) - 2j(j + 1)] (8) 

where the parameters a and 7 in (8) are determined by expressing this new Hamiltonian as 
an inverted SDI with the following choice 

1Q-j(j + 1) j(j + l)-3 
a - L, 2 H Tj W 

7 = 7^(C4-C 2 ) (9) 

Explicitly, they are 

20(j-2)(j + 3)(2j-3)(2j + 5) 
a j(j + l)(2j-l)(2j + l)(2 J + 3) 1 ] 

7 = 15[4j(j + l)-27] 

j(j + l)(2j-l)(2j + l)(2j + 3) 1 J 

Using the above parameters, one can show that the multipole pairing strength C' A are 

15(A-2)(A + 3)(A-4)(A + 5) 
A " j(j + l)(2j-l)(2j + l)(2j + 3) 1 ] 



The invertness can also explicitly be seen numerically. For j=21/2, the values of multipole- 
pairing strengths are 

[-0.001, 0, 0, -0.0093, -0.04, -0.11, -0.24, -0.46, -0.79, -1.27, -1.95] 

for A = 0,A = 2,A = 4, ,A = 20. 

Applying TDA to the inverted SDI, we can obtain the following energy expression: 

£ = E Wa + + l[L(L + 1) - nj(j + 1)] (13) 

A A 

Here the TDA energies,^ N A C' A are degenerate for all L values and thus serve as the 
intrinsic states from which the various rotational bands are built. The intrinsic excitations 
are the pairs with the energy C' A . The expression (13) with the numerical values of (12) 
clearly indicate that the pairs with the lowest energy does indeed have the highest spin, and 
hence the system's ground state must be due to the contribution of N pairs of nucleons, 
each coupled to the highest spin (A = 2j — 1). This is of course precisely the behavior 
demonstrated in a pair shell model calculation with the optimal pair approximation 3 . For 
the two pair case, the ground band is labeled as (2j — 1; 2j — 1) since it is formed by two 
2j — 1 pairs and the spin sequence of the band is 0, 2, 4, 6, while the first excited band 
is a (2j-3;2j-l)-band (one 2j — 3 pair and one 2j — 1 pair) with the 7-band spin sequence of 
2, 3, 4, 5, Finally, the /3-band comes from coupling two 2j — 3-pairs. 

In Table I, results for a system of four particles in a j = 21/2 orbit are presented. In 
particular, the above approximation is compared with the exact shell model calculations. 
In this approximation, the values of C' A are given by (12), referred to as Approx. I. This 
distinguishes it from the other choices of C' A mentioned below. In view of the simplicity of 
the model, the agreement is quite good. 

It is interesting to note that the above choice of C' A is close to that of fixing a and 7 in 
(8) which result in cancelling the a + 6A(A + 1) terms in the C\ for QQ interaction. This 
latter is given as 10 W(jjjj; 2 A) and can be expanded in A(A + 1): 

C A = a + 6A(A + 1) + cA 2 (A + l) 2 (14) 
5 



where 
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15(4j(i + 1) - 



j(i + l)(2j - l)(2j + l)(2j + 3) 
15 



j(j + 1)(2j-1)(2j + 1)(2j + 3) 
The a and 7 thus determined are 



80j(j + l) , , 

a (2 J -l)(2 J + l)(2j + 3) [b) 

7 = ^(^l)-!) (16) 

j(j + l)(2j-l)(2j + l)(2j + 3) 1 J 

and the new multipole pairing strengths are 

a = 15A2 ( A + 1 ) 2 fl7 ) 

A " j(j + l)(2j-l)(2j + l)(2j + 3) 1 ] 

which are numerically close to that given by (12) (See Table II). 

The above results imply that one can decompose the multipole pairing strengths in the 
QQ-interaction into the collective part a + 6A(A + 1) and the intrinsic part cA 2 (A + 1) 2 . The 
latter contributes to the intrinsic energies of TDA. The approiximation led by this choice of 
C' k is given as Approx.II in Table I. 

Besides the above choices of a and 7 which led to the inverted SDI-like interactions , 
an additional one is added in which the parameters are fixed so that C' A is as close to an 
inverted SDI as possible. A strict inverted SDI is defined as 

c "' = 27TT C 2-'-* (18) 

and the x 2- fit treatment would lead to the following expressions for a and 7: 

< (C A - Cl)A{A + i) > - < (C A - CD >< A(A + 1) > 
7 <A 2 (A + 1) 2 >-<A(A + 1)> 2 1 ' 
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a =< (C A -C a v )> - 7 [< A(A + 1) > -2j(j + 1)] (20) 



where < .... > represents the average taken as 



EaIq (2A + 1)(....) 
E(2A + 1) 

When these are applied to the QQ force, the resulting analytical expressions are complicated 
and not physically transparant. If one choose the following simple expressions (which will 
fit (19) and (20) well), 

8Q(j-3)(j + 4)(j-5)(j + 6) 
j(j + l)(2j-l)(2j + l)(2 J + 3) 1 } 

n, = 60(j-5)(j + 6) 

7 j(j + l)(2j-l)(2j + l)(2j + 3) 1 ] 

one can obtain the results as presented in Tables I and II (denoted as Inverted SDI). 

The comparisons we presented in Table I with the various choices of C' A suggest that the 
TDA with the inverted-SDI is an excellent caricature of the exact shell model calculations for 
the single-j case. Consequently, the intrinsic excitations induced by the QQ-interaction are 
described according to our simple model (13) as elementary with energies C' A , from which all 
the rotational bands are built upon, which are the multipole-pairing strengths in an inverted 
SDI. We therefore conclude from (8) that the "intrinsic" part of the QQ force is equivalent 
to an inverted SDI. With this view point, one can understand the origin of the optimum pair 
description of rotation. In physical terms, the short range delta force will favor the S-pairs, 
the QQ force favors high spins pairs whose field is generated by an inverted SDI. 

The present model given by Eq. (13) has two shortcomings. The first is the moment of 
inertias are identical for all bands and the second is that the parameter a is independent of 
the elementary excitations . Using a more rigorous treatment, (shown next) for the TDA for 
arbitrary interactions, the above results can be improved. In this treatment, which is akin 
to the generalized pair field method (GPFM) 4 , the results for the QQ interaction are nearly 
exact and more importantly, the moment of inertia can vary from band to band. 



We thus begin by considering the commutator of the auxiliary Hamiltonian 

H' = H-u;X^ (23) 

i<j 

and Ot in (4). H' can formally be expanded as 

[H',ft f ] = E' T rf + Yl $(12...iV;l , 2 / ...Ar';cj)(A t 1 ,A^....Aj v ,) L (24) 

l'2'...N' 

where u is chosen by a x 2_ fit with only the first term surviving. Again, the TDA is made 
valid by the proper choice of the new Hamiltonian! 

For a system of four particles moving in a single j-orbit, the u>, after x 2 -fitting, takes on 
the form 

= El^CClI + C L >)(2L[ + 1)(2L' 2 + l)X\ijL : :jjL,: L[L' 2 L) 
U 2E L/lL i 2 (2L[ + l)(2L' 2 + l)X^jjL 1]3 jL 2 ;L' 1 L' 2 L) { } 

where X is a 9j-symbol. Note that unlike a in (9), u is not a constant depends on L\L 2 L. 
This means that for the various bands, it will have different contributions to the elementary 
excitations. In addition, since it depends on L, uj can cover more than the term a + 7(A(A + 
1) — 2j(j + 1)) in (8) which is only limited to the rigid rotor case. 

For any general interaction C\, the u> in (25) is usually computed numerically. However, 
for QQ, we are able to find analytical solutions for the two pair case. 

Using C® Q = 10W(jjjj; 2A) for QQ, one can derive the following expressions: 

qq = IQgLi + 1)(2L 2 + l)F 2 + (2Li + l)(2L 2 + 1) E fc (-l) fc (2A: + l)C® Q F k 
" " l + 5 LlL2 -2(2L 1 + l)(2L 2 + l)X(jjL 1 ;jjL 2 ;L 1 L 2 L) 1 } 

where is given by 

F k = W(L 1 jL 1 j;jk)W(L 2 jL 2 j;jk)W(L 1 L 2 L 1 L 2 ;Lk) (27) 
+W(L l3 L 2 j;jk) 2 W(L 1 L 2 L 2 L 1 ;Lk). 

For the QQ-interaction, u QQ turns out to be simple in which the largest spin pairs are 
preferred to be chosen as optimum pairs. One can numerically show that both the second 
term in the numerator and the denominator, and the second term of Eq. (27) are all shown 
to be negligible. Thus, 
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i + s LlL . 



■(2L 1 + l)(2L 2 + l) 



(28) 



x W(L 1 jL 1 j-J2)W(L 2 jL 2 j-J2)W(L 1 L 2 L 1 L 2 ; L2) 



This expression will reproduce well the exact values for the low-lying bands. The detail 
comparisons with the exact results are given in Fig. 1. 

With this u, we can now cast the total energy into an analytical form: 



E = E + ^L(L + l)- X L 2 (L + lf 



(29) 



where E Q is the intrinsic energy given by 

Eo = cg? + c2?- x 

with the x i n the form 

15[3L! 2 -4j 2 -3][3L 2 2 -4P-3] 



— 4 — 4 — 2 — 2 2— 2~ 2 

Ll + L2 — L\ — -L 2 + g-^l "^2 



(30) 



j(j + l)(2j - l)(2j + l)(2j + 3)(2Lx - 1)(2L 2 - 1)(2L X + 3)(2L 2 + 3) 
and Lj = Li(Li + 1), j = yj(j + 1). The moments of inertia are found to be 

h 2 



(31) 
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X[2LY + 2LY - 1] 



(32) 



The term xL 2 [L + l) 2 which causes deviation from the rigid rotor expression is in fact 
negligible for the QQ-interaction. 

Since the choice of ui®® taken from a four-particle system is more accurate than those of 
a and 7 in (9), we now return to the rotational limit of the many particle system (13) and 
seek a better set of values them, extracted from Eq.(29). We thus set 



C QQ + C QQ + Aa _ 1 [Li(Li + 1} + L ^ L2 + 1}] = Eq 



(33) 



from which one can obtain a and 7: 



a = 



X 



— 4 — 4 10 — 2 — 2 



7 = 



2%~ X 



2Li + 2L 2 2 



(34) 



Since these expressions are extracted from the nearly exact solutions of the four particle 
problem, their ability to reproduce the inverted SDI-like behavior, as given by Eq.(8), is a 
confirmation of the physics so far discussed. Together with the former ones, the results of 
the moments of inertia and the various band heads, calculated according to this choice of 
parameters, are presented in Table I (listed as 2-Pair Fit). The values of C' A are listed in 
Table II. It is interesting to observe that these results are indeed close to the inverted SDI 
ones and the exact SM. Of course, it remains to be seen whether these expressions can be 
generalized to the system of an arbitrary number of pairs. 

In conclusion, we have presented evidence in this paper that the excitations induced by 
the QQ-interaction is equivalent to the inverted surface-delta interaction whose multipole- 
pairing component is dominant and attractive for the largest multipole. This thus provide an 
understanding of the origin of the recently introduced optimum pair description of rotation. 
Finally, it should be emphasized that the present approach possesses a predictive power for 
the excited bands using a simple Tamm-Dancoff approximation. 

Research reported in this paper is supported by the Chinese National Natural Science 
Foundation and the United States National Science Foundation. 
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Table I 



Inverted SDI vs EXACT SM 



Methods 


C'20 


a 


7 


7- E 


P-E 


Approx.I 


-1.95 


0.84 


5.4 x 10~ 3 


0.68 


1.36 


Approx.II 


-2.08 


0.91 


5.6 x 10~ 3 


0.70 


1.40 


Inverted SDI 


-1.53 


0.62 


4.3 x 10~ 3 


0.59 


1.18 


2-Pair Fit 


-1.49 


0.59 


4.2 x 10~ 3 


0.62 


1.12 


Exact SM 


-1.47 


0.58 


* 4.2 x 10~ 3 


0.62 


1.12 



** Note: In Exact SM and the 2- Pair Fit, the moments of inertia vary as follows: 

(a) fj =4.2 x 10~ 3 for g-band 

(b) g =3.3 x 10~ 3 for 7-band 

(c) g =2.5 x 10~ 3 for /?-band 
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Table II 



The Values of C' h 



A 


I 


II 


III 


IV 





-0.001 


-0 


-0.042 


-0.029 


2 





-4.E-4 


-0.034 


-0.020 


4 





-0.005 


-0.019 


-0.004 


6 


-0.009 


-0.021 


-0.019 


-0.014 


8 


-0.040 


-0.061 


-0.004 


-0.018 


10 


-0.110 


-0.142 


-0.003 


-0.006 


12 


-0.239 


-0.286 


-0.110 


-0.081 


14 


-0.456 


-0.519 


-0.267 


-0.232 


16 


-0.786 


-0.870 


-0.532 


-0.490 


18 


-1.273 


-1.376 


-0.939 


-0.891 


20 


-1.948 


-2.075 


-1.529 


-1.473 



Note. Column labels I, II, III and IV represent the following: 

I. Approximation I 
II. Approximation II 

III. Inverted SDI 

IV. Two Pair Fit 
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